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We review some recent developments in the theory of nonlinear von Neumann 
equations. We distinguish between the von Neumann equation (which can be 
nonlinear) and the Liouville equation (which should be linear). Explicit examples 
illustrate the technique of binary Darboux integration of nonlinear density matrix 
equations and special attention is payed to the problem of how to find physically 
nontrivial 'self-scattering' solutions. 



1 Why density matrices? 

"It is a common error always to view a mixed state as describing a system 
that is actually in one of a number of different possible pure states, with 
specified probabilities. While this 'ignorance interpretation' of the mixed state 
can indeed be a useful practical way to describe an ensemble of completely 
isolated systems, it entirely misses the deep and fundamental character of 
mixed states: If a system has any external correlations whatever, then its 
quantum state cannot be pure. Pure states are a rarity, enjoyed only by 
completely isolated systems. The states of externally correlated individual 
systems are fundamentally and irreducibly mixed. This has nothing to do 
with 'our ignorance'. It is a consequence of the existence of objective external 
correlation" EJ. 

The quotation from the Mermin paper can serve as a motto for what we 
are going to present below. A physical system that is described by a state 
vector at time t = will remain in a pure state for alH > if and only if it 
will never interact with anything. An interaction leads to correlations, and 
correlations mean non-product (entangled) states. A subsystem of a bigger 
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composite system which, as a whole, is in an entangled state, is no longer 
described by a state vector but by a non-pure reduced density matrix. The 
density matrix is an entirely quantum object whose ontological status is the 
same as this of the state vector it originated from. A theory that deals only 
with pure states deals either with the entire Universe or with objects that 
cannot interact and are therefore unobservable. 

The above remarks are especially relevant if one thinks of nonlinear gen- 
eralizations of quantum mechanics. 'Nonlinear quantum mechanics' is tra- 
ditionally associated with nonlinear Schrodinger equations. Such equations 
have the general form 

=F(V')|^). (1) 

Assuming we have a concrete physical system which, when isolated, evolves 
according to (j^) we immediately face the question of how to describe this 
system during or after an interaction. The problem has led to a great amount 
of confusion as to the real status of nonlinear generalizations of quantum me- 
chanics. Obviously, the resulting misunderstandings involve an interpretation 
of experiments testing linearity of quantum mechanics: A theory that does 
not tell us how to describe external correlations does not produce testable 
predictions and, hence, cannot be tested. 

The problem of subsystems versus (|l|) was, in our oppinion, solved by 
Polchinski in □. In short, his answer was the following: The equation describ- 
ing a composite system is 

= (i/3(ps)+i/r(Pr))|*) (2) 

where the index 's' stands for a subsystem and 'r' for the 'rest' or 'remaining 
systems'. The state vector represents the pure state of the composite 
'subsystem plus the rest' system, and the p's are the appropriate reduced 
density matrices. Let us note that, even if we assume, as we have done in (^), 
that the composite system is in a pure state, the formalism inevitably leads 
us to density-matrix-dependent nonlinear operators 
To include interactions one can write 

= (h,(Ps) + i?r(Pr) + i?int(«')) l^-) (3) 

with some interaction (linear or nonlinear) Hamiltonian operator Hinti'^)- 
The main feature of and (^) is the fact that for i/int(5') = (that is, 
when the interaction is over) the reduced density matrices satisfy the von 
Neumann-type equations 

iPs=[Hs{ps),Ps], (4) 

ip, = [H,{p,),p,]. (5) 
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The form of these von Neumann equations illustrates the locality of the 
Polchinski formulation: Subsystems that do not interact do not 'see' each 
other. It should be stressed that this property is typically lost when one 
considers different ways of extending the subsystem dynamics. An interest- 
ing discussion of the problem is given in these proceedings by Liicke a who 
shows that there may exist local extensions different from those discussed by 
Polchinski. However, the nonlinear Schrodinger equations they apply to are 
linearizable by a nonlinear gauge transformation □. The Polchinski-type ex- 
tensions bear some kind of universality but simultaneously are non-unique in 
the following sense: There exists an infinite number of inequivalent exten- 
sions which reduce to the same equation on product states-jV. more detailed 
discussion of this and related problems can be found in □'I'ErQu. 

From what we have written so far it follows that nonlinear von Neumann 
equations, whether we like it or not, will always appear when quantum non- 
linearity occurs. A somewhat more radical viewpoint is suggested by the non- 
uniqueness of the Polchinski-type extensions. To illustrate the point consider 
a simple spin- 1/2 nonlinear 'average energy' 

where CTz is the Pauli matrix. ThiS|-.type of nonlinearity was considered in 
experiments involving 2-level atoms cl. Now, how to write H{p) on the basis 
of H{i}i)l A natural guess is 

Hip) = (7) 



but 



Hip) ^ (8) 
f r p 

and an infinite number of similar fuctions would do as well. They all reduce 
to if p = However, if we start with some H{p) then the whole 

ambiguity disappears. This suggests that looking for a fundamental level of 
nonlinear quantum dynamics one should begin with von Neumann and not 
Schrodinger equations. 

Linear Schrodinger equation, as seen from a classical perspective, is sim- 
ply the equation of motion of an infinite-dimensional Hamiltonian system 
with average energy playing a role of Hamiltonian function. For a detailed 
presentatiorLof the formalism see the paper by Cirelli et al. published in these 
proceedings Linear von Neumann equation is also the equation of motion of 
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a classical Hamiltonian system with average energy in the role of the Hamil- 
tonian function. Geometrically, density matrices form an infinite-dimensional 
Poisson manifold with gUoo) Lie-Poisson bracket. This is described in the 
present volume by Bona 1121 and for this reason we will not spend much time 
on details of the Lie-Poisson formulation. 



2 Nonlinear von Neumann equation 

Let us denote by pa a matrix element of the density matrix p taken in some 
basis. The Lie-Poisson version of the linear von Neumann equation is 

tpa^{pa,{H)} (9) 

where (H) = Tr Hp is the Hamiltonian function. Nonlinear quantum mechan- 
ics d la Bona and Jordan W3 is based on the same Lie-Poisson structure but 
with (H) replaced by a nonlinear Hamiltonian function iJgon- The dynamics 
given by or its nonlinear generalization 

ipa ^ {Pa,Hgcn} (10) 

defines a Hamiltonian flow on the Poisson manifold of states. Having such a 
flow one can consider a classical probability density and its associated 
Liouville equation 

iw = {w,Hgcn}- (11) 

The classical probability density has a clear physical meaning: It describes 
a classical lack of knowledge about the state of a quantum source. In any 
experiment one faces this type of clasical uncertainty and all experimental 
averages one measures in a lab are of the form 

(A)exp = J dpw{p)TTAp (12) 

with J dp meaning an integration over the parameters defining initial con- 
ditions for the quantum dynamics. The Liouville equation ( pT| ) is linear in 
w independently of whether the von Neumann equation ( p^ is linear in pa- 
Moreover, the probability density w is directly accessible to the experimen- 
talist and reflects the classical configuration of the experimental setting. For 
this reason it is very important that the Liouville equation is linear. Formula 
12|) shows that it is practical to introduce the object 

Pcxp = / dpw{p)p (13) 
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satisfying 




Tr Ap< 



'cxp 



(14) 



and playing a role of a semiclassical density matrix. The 'common error', men- 
tioned by Mermin in the quotation we have started with, is the belief that all 
density matrices one encounters in quantum mechanics have such a semiclas- 
sical origin. The 'truely quantum' density matrices one obtains by reduction 
of entangled states to subsystems can be written in different bases of 'pure 
states' and all such bases are regarded as physically equivalent. Putting it 
differently, no particular decomposition of such a p into a convex combination 
of projectors is physically special. |_This is one of the important impossibility 
principles of quantum mechanics til. On the other hand, the decomposition 
defined by w is not only very special but, actually, is even uniquely given by 
the form of the experiment. An experimentalist can arbitrarily tamper with 
w but different convex combinations forming a concrete p are definitely out 
of his reach. 

One can imagine also physical situations which are somehow in-between 
the cases we have just described. For example, consider an entangled pair of 
particles and an experiment where we have a shutter which is opened whenever 
a particle labeled '1', say, is measured and a concrete result (say, spin 'up' 
or 'down' in a 'z' direction) is found. Each time the shutter is opened the 
particle labeled '2' leaves a box and, hence, the box is a source of particles 
in a concrete mixed state which depends on the entangled state of the pair. 
The mixed state, as depending on the macroscopic and clasically controlled 
actions undertaken by the experimentalist, is no longer 'fully quantum'. The 
resulting mixture is of a pexp type and there is no reason for the corresponding 
dynamics of the density matrix to be nonlinear. It seems reasonable to assume 
that nonlinear quantum dynamics of mixed states can occur only in cases 
where the very form of the 'pure-state' decomposition is in principle out of 
control. Such an impossibility principle seems to eliminate all the problems 



The distinction between the von Neumann equation ( pi3| ) and the Liouville 
equation ( pT| ) is therefore essential. One should not use the misleading term 
the 'Liouville-von Neumann equation' which suggests that the von Neumann 
equation is simply a quantized version of the linear Liouville equation and, 
accordingly, must be linear as well. 

Historically, the linearity of the standard von Neumann equation seems to 
have its roots in the linearity of the Schrodinger equation. The two equations. 



analyzed by Polchinski in 



(15) 
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combined with p ~ |V'fc)(V-'fc| imply 

ip^[H,p]. (16) 

Since the equation is linear, one can consider more general convex combina- 
tions 

P = Y,M^k){M (17) 

k 

which again satisfy (p^. This simple argument looks so natural that one 
may have not even noticed the additional assumption we have smuggled in. 
Indeed, we have started with (|l^) which was valid for pure states, that is those 
satisfying p^ = p, and by linearity extended the argument to p ^ p^ . But a 
pure p has eigenvalues and 1. Therefore any function / satisfying /(O) = 
and /(I) = 1 will satisfy f{p) = p if p is pure! It follows that the linearity of 
the Schrodinger equation implies at most an equation of the form 

ip=[HJ{p)] (18) 

with f{p)] — [H,p\ for p^ — p. In the next sections we shall devote much 
attention to the so-called Euler-Arnold-von Neumann equation obtained if 

f{x)=x\ 

Now, what is the relation between (|l^) and (p^? The answer is rather 

arprising: ([ij) is an example of ([l0| ) with an appropriate choice of -ffgen 
. To see how this works assume that for x £ [0,1] the function / has a 
convergent Taylor expansion 

oo 



fc=0 



and define 



i7gen(p) = Tr(/(p)i/) (20) 

(in the context of 111 it is important that Hgcn{p) is 1-homogeneous; this 
complication is left out here). Then ( [lO| ) is equivalent to 

ip^[H{p),p] (21) 

where, the effective nonlinear Hamiltonian is defined by the functional deriva- 
tive Il3 

H{p) = yH^M = T.fkT.(P- cil)''-'-"H{p - air (22) 

^ k=l n=0 
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and 

ip^[H{p),p]^[H,f{p)]. (23) 

The form implies that C„ — Tr (p") is time- independent for any natural 
n (such Cn are Casimir functions for {•,•}) from which iL. follows that the 
spectrum of a solution p of ( |l8|) has to be time-invariant E3. As we can see, 
although the dynamics of p is nonlinear for non-pure p, we are nevertheless 
still quite close to linear quantum mechanics. Looking at the explicit solutions 
for f{x) = we will see that we are, in fact, surprisingly close to the linear 
dynamics but with very specific and subtle new nonlinear effects at hand. 

It is quite clear there are fundamental reasons for investigating nonlin- 
ear density matrix equations. But one does not have to believe in nonlinear 
quantum mechanics to investigate nonlinear evolutions of mixed states. It is 
known that nonlinear voii Neumann equations arise very naturally in various 
mean-field approaches lij where they are solved either numerically or by ap- 
proximate methods. The program we are developing aims at finding exact 
techniques of solving nonlinear von Neumann equations. 

3 Lax representation of von Neumann equations and binary 
Darboux covariance 



The class of equations we have under some control is of the form 




n n 

ip = [H{p),p] = = Y.[A-\pA^p] (24) 

where A is a time-independent self-adjoint operator. For n = 1 and A — H 
one finds 

H{p) =Hp + pH (25) 

and the equation becomes the Euler-Arnold-von Neumann equation 

ip=[Hp + pH,p] = [H,p^]. (26) 

For all n the von Neumann equations involve Hamiltonians H{p) which are 
linear in p, and therefore the nonlinearities of the equations are quadratic. 
This property is important as it seems that only the quadratic nonlinearities 
can be treated by the binary Darboux transformation. In fact, what makes 
the class so interesting is the knowledge of its binary-Darboux-covariant Lax 
representation 

z^\ip) = ip-pAM, (27) 
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n 

= (28) 



fe=o 



where z^, ^ e C. The necessary condition for the solution \ip) of (P7|)-(P8|) 
to exist is that the relation (|2j) between A and p holds true. The idea of 
the Darboux transformation is, given a density matrix p and a solution 



of (|27|)-(p8[), to find a new solution |</5[1]) and a new p[\\ which are again 
related by the pair (|27|)-(|2^) with A unchanged. The fact that is a 

solution means that the ncsessary condition for its existence is satisfied, and 
this implies that p[l] and A are again related by (|2^). One may think of p and 
p[V\ as 'ground' and 'first excited' states, and the Darboux transformation is 
a 'creation operator'. 

The so-called binary Darboux trpisformatinn is constructed as follows (for 
more details and generalizations cf. lla^EJ^IlZH13lij). Take some linear operators 
V and J, a parameter s, and three complex numbers p, v, and A. Now 
consider three linear equations 

idsW)^{V-pJM, (29) 
~idM^{x\{V-vJ), (30) 
-id,{',p\^ {■<P\{V -\J). (31) 

It is important that ( pO[ ) and ( |3l] ) involve bras and (|2^) involves a ket. Define 

P = MM (32) 

V[l] = V + {p-u)[P,J] (33) 
(^[1]H(V'|(1-^P) (34) 

(the fact that we transform here (^| and not \ip) or (x| is not essential; one 
can construct similar transformations of any of these solutions with the help 
of the remaining two). A straightforward calculation then shows that 

-zd,{m\ = {m\{v[i]-\j). (35) 

To apply the above technique to our von Neumann equations we take three 
pairs: (p7|)-(p8|), which we already have, and two more with parameters A, za, 

za(V'I = (V'I(P-AA), (36) 

n 

= (V^l ( ^ A"-^^^ - AA"+i) , (37) 

k=0 
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zAx\ = {x\{p~^A), (38) 

-m = (xi ( E ^""V^' - ^^"+') ■ (39) 



fc=0 



In what follows the and J^-pairs (27)-(|2^) and (|3q)-( |39|) w ill be used to 
define the binary transformation of the conjugated A-pair (p6[) - (pTj) . 
Defining P as above and 

p[l]= p+{^i-v)[P,A] (40) 

we find 

z>^{m\ = {m\m-^A) (4i) 

n 

-z{m\ = (^[1]|(EA"-V[1]^'=-AA"+i). (42) 

fe=0 

If all the objects necessary for the construction of exist, then the nec- 

essary condition for its existence must be fulfilled, and this means 

n n 

Having one solution p we have managed to produce another solution p[l]. In 
this respect the Darboux transformation is a really wonderful device. How- 
ever, from time to time surprises can occur. It is clear that if all the assump- 
tions we have made are satisfied then p[l] must be a solution. But the general 
theorem does not guarantee that the solution is nontrivial! In fact p = A, 
p = 1, or even p — are also solutions of (p^), so can we guarantee that such 
pathological cases are excluded and p[l] is physically interesting? 

The answer depends on what is meant by 'interesting'. In general, each 
Darboux transformation has an inverse. In classical problems, such as the 
Korteweg-de Vries equation, it is typical to start with a trivial solution u = 
and u[l] is already a soliton, definitely a highly nontrivial solution E3. This 
means there exists a transformation that maps a soliton into 0. 

We will now show that this type of pathology is excluded if a binary trans- 
formation of the type we use is considered. What is not excluded, however, 
are the cases when, say, p[l] — p or p[l] = T pT^^ , with T a time-independent 
unitary transformation. Such a possibility leads to nontrivial practical com- 
plications. The next theorem is of fundamental importance and its proof is so 
elementary that we can give it here E£l. Consider the following three general 
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Lax pairs 

^{P~ l^A)\^), (44) 

^\^)^{V{p)^^lJ)\^), (45) 

zAx\ = {x\{p->^A), (46) 

-^{x\^{x\{V{p)-vj), (47) 

ZAiV'l = (V-Kp-AA), (48) 

= (^l(y(p)-Aj). (49) 

The only assumption we make about p, J, A, and V{p) is the covariance of 
(^8|)-(^9|) under the binary Darboux transformation constructed with the help 
of (xl and \ip). 

Theorem 1. Under the above assumptions the binary Darboux transforma- 
tion p 1-^ p[l] is a similarity transformation, p[l] — TpT^^, where 

r=l + ^P = e^'"^ (50) 

Proof: By definition p[l\ = p + [p - v)[P, A]. Eqs. (H), (||) imply 

z^P ={p- pA)P (51) 
z^P = P{p - vA) (52) 

from which it follows that 

P{p^pA)P={p-pA)P (53) 
P{p-vA)P = P{p-vA). (54) 

Multiplying (js^) by ( ^) by /i, and subtracting the resulting equations we 
get 

[P, A] = "^—itppP --pP+-Pp. (55) 
pi/ p V 

Inserting this expression into the definition of we obtain 

p[l]^(l + ^p)p(l + ^p). (56) 

Q.E.D. 

The theorem explains why p\\\ — etc. is excluded if p is a density matrix: 
Spectra of p and p\\\ must be identical. This property is characteristic of all 
equations that are compatibility conditions for binary-Darboux-covariant Lax 
pairs. The result will, in general, not hold if one considers a binary-Darboux- 
covariant zero-curvature pair of the type used for example in Ej. 
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The next important result, whose proof follows from a straightforward 
calculation, is the covariance under 'spectnmi shifting and rescaling'. 

Theorem 2. Assume [X, A] = [X, p] = 0, F G R, and p is a solution of (|2|). 
Then 

px{t) = e-*("+i)^^"*(p(t) + x)e*("+i)^^"* (57) 
PY{t) = YpiYt) (58) 

also satisfy (U). 



4 Two strategies leading to 'interesting' solutions 

The problem with the Darboux transformation is that in order to find a 
solution p[l\ we must already know somehow another solution p to start with. 
As we have seen above, there are some obvious solutions such as 0, but by 
Theorem 1 they will not lead to anything nontrivial. Still, the case is not 
completely hopeless. 

Ji^.l The first strategy 

The strategy works for the Euler- Arnold- von Neumann equation El 

ip^[H,p\ (59) 

In all the examples discussed below we assume v = p, and \ip) = \x) since this 
guarantees that hermiticity is conserved by the transformation. A pure state 
p = p^ \s simultaneously a solution of ( p^ ) and of the linear von Neumann 
equation with the same H . Therefore, pure state solutions of the linear equa- 
tion form a nontrivial subset of solutions of the nonlinear one. Unfortunately, 
by Theorem 1 we will have p[l\ — p[\Y and although such cannot be 
claimed 'trivial' they are nevertheless quite 'uninteresting'. Interesting solu- 
tions are obtained if one starts with p satisfying [iJ, p^ — ap] = 0, for some 
a G R, but sucb_that the operator := p^ — ap ^ Q is not a constant times 
1. In this case 113 

p{t) = e-"'^V(0)e"'-^* (60) 

and 

p[m = e-^'^^*(p(0) + (M-Ai)i^a(t)-^e-^^»*[|^(0))(^(0)|,i7]et^"*)e-^*, 
=: e-''^«Vint[l](i)e"^* (61) 



11 



where is a complex parameter of the Darboux transformation, |v(0)) a 
solution of the Lax pair at t = 0, and 



F,(t) = (^(0)|exp(z^A,t)|<^(0)) 



Ml 

After this first step we know were to look, but we still have to find an ap- 
propriate p and \ip). It turns out we need three more tricks which are best 
illustrated by the following example (see the Appendix). 
Consider the Hamiltonian 



if = 1 (62) 




and take fi = i (for real fi the binary transformation is trivial). We have to 
find an appropriate initial condition p{0). The first trick we have mentioned 
is to begin with a solution which is neither normalized nor positive. Such a 
non-density-matrix solution will be denoted by ^ instead of p. We can always 
make it positive and normalized by using the transformations of Theorem 2. 
We take ^{t) = e-'"*^{0)e'"* with 



(63) 



The first two matrix elements on the diagonal are the solutions of the equation 
— X = 1/A (this is the second trick) and for this reason 





which obviously commutes with H (1/4 is not essential here; this could be 
basically any nonzero number). Therefore [H,£^{t)'^] = [H,^{t)] even though 
^(i)^ 7^ C(^) ^^^d [H,£^{t)] 7^ 0. It may happen that appropriate solutions 
of the quadratic equation x^ ~ x — xq do not exist. But to make the trick 
work it is sufficient to find solutions of x'^ — ax — xq with some a and then 
we will have Aa instead of Ai. The third trick is to choose ^(0) and \(p) in 
such a way that the solution of the Lax pair is not an eigenstate of A^ (since 
the contributions from Fa{t), exp(— ^Aai), and exp{j-Aat) would cancel one 
another in (|6l| ) , and the internal part of (|6ll ) would become time-independent 
— that is exactly what we want to avoid). 
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The eigenvalues of ^(0) — iH are z± = (1 ± 'i\/2) /2 and z_ has degeneracy 
2. The two eigenvectors corresponding to Z- are orthogonal: 





Taking an arbitrary linear combination of them, say, 



|^(0)) = i=(|^i) + |¥'2)) 



(65) 



(66) 



we obtain an eigenvector all the three components of which are non- vanishing. 
This implies that the unitary transformation T occuring in Theorem 1 will not 
be block-diagonal. The fact that the two orthogonal eigenvectors correspond 
to the same eigenvalue of ^ — /ii? and to different eigenvalues of allows 
us to construct a nonlinear dynamics involving the entire 3-dimensional space 
and automatically guarantees that \(p) is not an eigenstate of Ao. We find 



with 



I l+\/2 _ V2 

2 1+e' 



\ 2\/2cosh(t/2) 



l-\/2 



2cosh(t/2) 



2i/2cosh(t/2) 

72 1 
1 + e' 2cosh(t/2) 

2 



(67) 



(68) 



^(t) is not yet a density matrix since it has a negative^ eigem ^-alue (1 - V2)/2 
and its trace is 3/2. Now it is time to use transformations from Theorem 2. 
Set X = "^"^ 1 and Y = \/2/3. A combination of the two transformations 
leads to the final solution 



where 



Pint(i) = 



PxY[im = e-^3«Vint(t)e' 



/ V2- 



V2 



l+eN/2t/3 







V2 



■Ht 



2\/2cosh[t/(3y2)] 
1 



l+e^*/3 2cosh[i/(3\/^)] 

, ^ 1 vl / 

\ 2\/2cosh[t/(3v^)] 2cosh[t/(3\/2)] 2 / 



(69) 



(70) 



The density matrix pxy[1](0 has eigenvalues 2/3, 1/3, and 0. For \t\ — > oo 
one gets the dynamics asymptotically linear but with different asymptotics 
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Figure 1. Projection of the average spin on the x-y plane for times < t < 10 (in dimen- 
sionless units). The amphtude of the oscillation decreases with time. 



for negative and positive times, and a kind of 'self-scattering' (or 'phase tran- 
sition' as we called it in E3) around t = 0. To have some qualitative feel of 
what happens consider averages of spin-1 matrices 



The figures illustrate the effect. Figures 1 and 2 show the evolution of the 
x-y projection of (J). The magnitude of the oscillation goes very quickly 
to as \t\ grows. Figure 3 shows the dynamics of the z component of (J). 
The process of 'self-scattering' maps one asymptotically linear solution into 
another. A discussion of more realistic systems, including a one-dimensional 
harmonic oscillator, can be found in Ej. 

4.2 The second strategy 

The first strategy leads to solutions involving a finite number of states and 
cannot be directly applied to n > 1 equations from the Darboux covariant 
family (H): The non-abelian shape of, say, the n = 2 equation 




(71) 
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3x10 
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2x10 
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^arxio 


23 






-•^xio"^-^x]?r^ 




5xl0"^*Lxl0'i^.5xl0"" 



Figure 2. Seeds of destruction of the past asymptotic state: The same dynamics as in Fig. 1 
but for times -230 < t < -220. The amplitude is 10^^ times smaller than at Fig. 1. The 
amplitude of the oscillation increases with time. 




Figure 3. The 'self-scattering'; Average (Jz) as a function of time (solid). The dynamics 
is asymptotically linear but with different asymptotics for t — > +00 (dashed) and t — » —00 
(dotted). The self-scattering phase shift is clearly visible. 
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requires different tricks. The solutions discussed in the previous subsection 
started with non-stationary p's, since those satisfying [p, A] = lead to the 
projector P commuting with both A and p, and the binary transformation is 
trivial. Still, there exists another class of stationary solutions of (p^, obtained 
if Ap = —pA. Now the projector P will not, in general, commute with p and 
A, and the binary transformation may be nontrivial. 

Assume [P, p] ^ 0. For n = 2m one has X]fe=o(^-'^)'^ ^ ^ ^^'^ 

»|^) - (a>-/M"+i)|^) = z^A^l^) (72) 
The formal solution of the latter equation is 

|^(i))=e--^"*|^(0)). (73) 
For odd n we get J2k=oi^^)'' = ^ ^^'^ 

m = -/iA"+V> (74) 

implying 

|^(i))=e^M"+^*|^(0)). (75) 

To have an 'interesting' p[l] we must make sure that the solution of the Lax 
pair is not an eigenstate of A^"^ (for n = 2m) or yl^Cm+i) (£pj. ^ _ 2m + 1); 
otherwise we have a problem similar to this with the eigenstates of A^. Let us 
note that p which anticommutes with A does commute with A^ and therefore 
also with the generators of the time evolution given by (|7^) and (|7^) , and we 
can look for an eigenstate of p — pA at i = 0. 

For simplicity we shall check the trick again on the Euler- Arnold- von Neu- 
mann equation. Let ak and Cfc be Dirac-a and Pauli matrices, respectively. 
Consider the Hamiltonian 

J? = ai (g) l2x2 + l4x4 fXi 0-1, (76) 

and take a non-density-matrix stationary solution 

^ = a2 "JS 0-2 + Q!3 (8) (73 (77) 

satisfying £_H = —H^. Set p = i and take 

(^1 = (i, 0,-1,0, -z, 0,1,0) (78) 
satisfying {p — iH)\(f) ~ 0. Then 
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The eigenvalues of ^ and •CllJlO ^^'^ and ±2. To produce the density matrix 
we shift the spectrum by A > 2 and rescale to get the unit trace. This will 
be again a 'self-scattering' solution but its explicit form will not be shown 
here. The second strategy has the advantage of being applicable to generi- 
cally infinite-dimensional problems. A nontrivial question (requiring still more 
tricks) is how to produce trace-class solutions if the Hilbert space is not finite- 
dimensional. 
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Appendix 

Generic solutions of the equation of motion ( p9|) with H given by (^2|) can be 
obtained by a straightforward integration. Change the basis so that H is of 
the form 

/ 11 \ 
H = -[ -fi (79) 
\0 OA/ 

Then the diagonal elements of p are constants of motion. Write down the 
equation of motion for the remaining matrix elements. One obtains 

ipi,2 = 2^[(pi4 -I- P2,2)Pl,2 + Pl,3P2,3], 
*Pl.3 = - A)[(pi,l -I- P3,3)Pl,3 + Pl,2P2,3], 

j/02,3 = -(m + •^)[(P2,2 + P3,3)P2,3 + PTJPl,3]- (80) 

Let W — |pi.2p- It is not difficult to show that it satisfies the equation 

W = aW^ +bW + c (81) 
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with a, b, and c constants. Equation (81) is well known O. Its solutions are 
doubly-periodic functions. The result is of the form 



Wit) ^ p-'^sn^{a{t-to),k)+j (82) 

with sn the Jacobi elliptic function, and with fc, a, /3, 7, and to constants. In 
the limit A: = 1 one of the periods becomes infinite. The result ( |69| ) obtained 
by the binary Darboux transformation corresponds precisely to this k — 1 so- 
lution. It is at the moment not clear to us what kind of a seed solution (if any) 
can lead, via the binary Darboux transformation, to the k ^ 1 class. An al- 
gebraic characterization of such additional seed solutions would be important 
from the perspective of more general cases, especially the infinite-dimensional 
ones. 
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